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The gap equation for fermions in a version of thermal QED 
in three dimensions is studied numerically in the Schwinger- 
Dyson formalism. The interest in this theory has been re¬ 
cently revived since it has been proposed as a model of 
high-temperature superconductors. We include wave-function 
renormalization in our equations and use a non-bare vertex. 
We then have to solve a system of two integral equations 
by a relaxation algorithm. Fermion and photon self-energies 
varying independently with energy and momentum are used, 
which should produce more accurate results than in the pre¬ 
vious literature. The behaviour of the theory with increasing 
temperature and number of fermion flavours is then carefully 
analyzed. 


tion effects and a non-bare vertex. 

This complicates the study considerably and consti¬ 
tutes a non-trivial step forward, since instead of having 
only one gap equation, a system of two integral equa¬ 
tions has to be solved. We continue to neglect the imag¬ 
inary parts of the photon polarization functions and of 
the fermion self-energy for simplicity, as other studies 
have done so far & & Their inclusion, even though 
in principle necessary, would double the number of cou¬ 
pled equations to be solved, which would render the nu¬ 
merical study too complicated for our algorithm and the 
computer power at hand. 


I. THE PROBLEM 

The study of QED in three dimensions during the past 
years has revealed very interesting physics in this the¬ 
ory §- §• At finite temperature, a version of QED , 
i.e. t 3 — QED , has been proposed as a model of high- 
temperature superconductivity [Q. The behaviour of 
the theory with increasing temperature T and number 
of fermion flavours Nf was recently studied numerically 
in the Schwinger-Dyson formalism in a bare-vertex ap¬ 
proximation This approximation corresponds to a 
particular truncation of the Schwinger-Dyson hierarchy, 
and its effects, along with wave-function renormalization 
effects, could be important |jj, especially in connection 
with the behaviour of the theory when the number of 
fermion flavours is varied. This has led to several in¬ 
teresting studies at zero-temperature trying to include 
these effects 0- in order to solve the controversy on 
whether there is a critical number of fermion flavours 
beyond which there is no mass generation. At finite tem¬ 
perature, similar studies have so far used several approx¬ 
imations that we will relax in this study. 

The first purpose of this work is to provide an accurate 
value for the ratio r = 2M(0, 0 )/ksT c Q, where M(0,0) 
is the fermion mass function at zero momentum and en¬ 
ergy and T c is the critical temperature beyond which 
there is no mass generation. Apart from being theoreti¬ 
cally very interesting, this quantity is of direct physical 
relevance since it can be compared to corresponding val¬ 
ues measured in certain high-temperature superconduc¬ 
tors. The second purpose is to provide a reliable phase 
diagram of the theory with respect to T and Nf. We 
continue to consider momentum- and energy-dependent 
self-energies, but we include wave-function renormaliza¬ 


II. THE EQUATIONS 


We will employ the Schwinger-Dyson formalism in or¬ 
der to gain information on the momentum-dependent 
fermion self-energy and the non-perturbative physics be¬ 
hind it. The Schwinger-Dyson equation for the fermion 
self-energy is given by 
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where q = p — k, e is the dimensionful gauge cou¬ 
pling of the theory which we will take to be constant 
throughout this study, A ^ is the photon propagator with 
p, v = 0,1, 2, T u is the full photon-fermion vertex, 7^ is a 
four-dimensional representation of the 7-matrices, So is 
the bare fermion propagator, and the finite-temperature 
fermion propagator in the real-time formalism is given by 
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where /3 = 1/ksT, A{p) is the wave-function renormal¬ 
ization function, 5 is the usual Dirac function and we have 
made a rotation to Euclidean space. Note that we avoid 
the matrix form that the propagator has in this formal¬ 
ism, since the Schwinger-Dyson equation that we have 
written down involves only a one-loop diagram directly, 
so complications due to the field-doubling problem do 
not arise jf3|. However, it has to be noted that a more 
careful treatment involving the imaginary parts of the 
self-energies would involve the full matrix propagators 
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B . The resulting equations would in that case be unfor¬ 
tunately too complicated to be solved even numerically, 
without serious truncations. Moreover, due to the broken 
Lorenz invariance at finite temperature, the wave func¬ 
tion renormalization could in principle affect differently 
the po and \p\ propagator components, i.e. we should re¬ 
place (1 + A(p))j> by ((l + A(p))'y°+a)p 0 + (l + B(p))'y*pi 
with i = 1,2. For simplicity we will restrict ourselves 
to situations where a = 0 , which correspond to a zero 
chemical potential, and we will work in the approxima¬ 
tion where A(p) = B(p) also for non-zero temperatures 
as done in similar studies 

For the vertex T u {p, q) we use the ansatz T v {p, q) = 
(1 + A(q)) 7 ", where A(q) is the same wave-function 
renormalization function appearing in the fermion prop¬ 
agator. It is unfortunately not symmetric in the vertex 
momenta, but its use simplifies the numerical algorithm 
considerably. Even though this vertex does not satisfy 
a priori the Ward-Takahashi identities, it is expected 
to incorporate the basic qualitative features of a non- 
perturbative vertex at zero temperature when used in a 
Schwinger-Dyson context |10. It has also been used in a 
finite-temperature case jl 2 | 7 "supported by the qualitative 
agreement of the results of this ansatz with the ones ob¬ 
tained by more elaborate treatments 0 . Furthermore, in 
the problem of the present study it gives results similar 
to a symmetrized vertex 

Moreover, the photon propagator in the Landau gauge 
is given by @ 
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with i,j = 1 , 2 , and where we neglect its temperature- 
dependent delta-function part since it is expected to give 
a vanishingly small contribution ||, [jllj] , JlTj . The longi¬ 
tudinal and transverse photon polarization functions II P 
and II p are given explicitly in 0 and taken from [ pT[ , 
where they are calculated in a massless-fermion approx¬ 
imation and where wave-function renormalization and 
vertex effects cancel only for specific vertex choices 0 . 
One should in principle couple the expressions for II pp in 
our system of integral equations in a self-consistent man¬ 
ner, but unfortunately this would render our numerical 
algorithm too complicated. 

Identifying the parts of this equation with the same 
spinor structure, we reduce the problem to that of a sys¬ 
tem of two three-dimensional integral equations involving 
two functions varying independently with po and \p \. The 
equations take the following form: 
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where a = e 2 Nf, and it is more convenient to work with 
the mass function M(p 0 , \p\) = E (p 0 , |p|)/(l+A(p 0 , |p|)). 
We also sum over the photon polarizations P = L,T and 
over the two roots of the delta function by introducing 
e = 1,-1. The quantity E is approximated by the re¬ 
lation E 2 ss |q| 2 + M 2 (0,0) S, where use of the delta- 


function property 8{ax) = 8{x)/\a\ has been made. 
Furthermore, the functions Q and P are given by 
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One can easily check that Q + P = —2 {pk){kq)/k 2 , which 
would reproduce the result of |l2j if one takes II P (A:) = 
II p (fc) = II(fc) and switches to imaginary-time formalism. 

After inspecting the equations, we note that the vertex 
ansatz we chose makes the integral giving the function 
A{po-i |p|) depend only on the function M and indepen¬ 
dent of A. On the other hand, the equation for M(p 0 , |p|) 
has to be solved self-consistently, and it actually always 
accepts, apart from the solutions we will seek, the trivial 
solution as well. An analytical study of such a system 
would not be possible without severe approximations, so 
we proceed to the numerical solution of the equations 
given above. 


III. NUMERICAL RESULTS 
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The numerical method used to solve the system of 
equations is the same as the one presented in [ 0 . The 























physical quantities of interest do not vary substantially 
with our ultra-violet (UV) cut-off, since at the effective 
UV cut-off a of the theory they are already negligibly 
small [^b). The solution at T = 0 and Nf = 2 for the 
functions E(p 0 , |p|) and — A(p 0 , |p|) is given in Figs. 1 
and 2 respectively. The general variation of these func¬ 
tions with momentum does not change with increasing 
temperature or varying of Nf, even though the overall 
scale of E drops fast with Nf. The function E (p 0 , |p|) 
falls as expected with increasing momentum, and is of the 
same form as the mass function M(po, |p|). The function 
A(po, |p|) is always in the range between -1 and 0 as re¬ 
quired |lg| ], it is tending to zero for increasing momenta, 
and it is of the same form and magnitude as the approx¬ 
imate form used in 0 - 

For a given number of fermion flavours Nf, when the 
temperature exceeds some critical value there is no solu¬ 
tion for the fermion mass function but the trivial one. 



FIG. 1. The fermion self-energy E(po, |p|) at zero temper¬ 
ature and for Nf = 2, A uvIcy. = 0.1, as a function of energy 
and momentum in logarithmic scale. All quantities are scaled 
by a. 



FIG. 2. The opposite of the wave-function renormaliza¬ 
tion — A(po, |p|) at zero temperature and for Nf = 2, 
A uv/ot = 0.1. 
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FIG. 3. The fermion mass function at zero momentum and 
zero temperature, scaled by a and on a logarithmic scale, with 
respect to Nf for a ratio A uv/a = 0.1. We fit our results with 
the curve M(0,0)/a = e ~ 1ASN f /13.25. Values of Nf larger 
than 3.35 are not considered, because then the self-energy 
falls below the IR-cut-off. 
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FIG. 4. The phase diagram of the theory with respect to 
the temperature scaled by a and in logarithmic scale and the 
number of fermion flavours, for a ratio A uv/a = 0.1. We fit 
the data points with the curve ksT/a = e ~ 1A8N f /70. This 
curve should not be extrapolated for Nf > 3.35. 

The value of the ratio r is concentrated approximately 
around r ss 10.6, which is comparable to the values ob¬ 
tained in H which neglected A(po,|p|)- The ratio r 
found is also comparable to typical values obtained in 
Ref. [^2), which includes wave-function renormalization 
effects but uses several approximations which we were 
able to by-pass in this study. We confirm therefore that 
adding these effects does not influence the behaviour of 
the theory in a significant way. We have to note more¬ 
over that our r values are somewhat larger than the value 
r w 8 measured for some high-temperature superconduc¬ 
tors Jll|. However, we could be overestimating this ratio 
because of a possibly poor convergence of the algorithm 
for temperatures close to the critical one. 

At zero-temperature, this theory is known to ex¬ 
hibit also an interesting behaviour with the number of 
fermions Nf. In Fig. 3 we plot the zero-momentum 
and zero-temperature fermion mass function with re¬ 
spect to Nf. We fit the data with the exponential curve 
e -1.48JV// 13 25 . 

At Nf ~ 3.35, the mass function is still roughly four 


No mass generation 
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times larger than the cut-off. When Nf > 3.35, our algo¬ 
rithm does not converge and the mass function tends fast 
below the IR cut-off. This behaviour could indicate that 
Nf « 3.35 is some critical point beyond which dynamical 
mass generation is impossible. 

The value of Nf we find is remarkably close to the 
one quoted in the numerical study of |J, and it is also 
close to our previous result (|]. A similar study [0 , which 
includes a calculation of the fermion-field anomalous di¬ 
mension to second-order in 1 /Nf, predicts a critical value 
Nf 3.28, which is only slightly larger than the theoret¬ 
ical prediction neglecting wave-function renormalization 
which gives Nf = p- ~ 3.24 0, and quite close to the 
value we find numerically. 

In Fig. 4 we plot the phase diagram of the theory with 
respect to Nf and fcgT. It separates two regions of the 
parameter space which either allow or do not allow dy¬ 
namical mass generation. The choice of an exponential 
fitting curve was only made to describe “phenomenolog¬ 
ically” the general tendency of the data and to provide 
a measure for a r-ratio independently of Nf , and is rem¬ 
iniscent of the results in Ref. || but with a somewhat 
steeper slope. However, there are also studies that pre¬ 
dict a non-analytic behaviour of £ for Nf near its critical 
value 0. Lack of convergence of the algorithm and fall 
of the mass function below the IR cut-off however do not 
allow us to test the precise behaviour of the theory near 
Nf ss 3.35, and it is also clear that the fitting curves 
should not be extrapolated for Nf larger than this value. 

IV. DISCUSSION 

We were able to solve a system of two coupled inte¬ 
gral equations for the fermion mass function and wave- 
function renormalization for a finite-temperature version 
of three-dimensional QED , by applying a numerical re¬ 
laxation technique. One main result is a r-ratio of about 
10.6, which is close to previous numerical studies, con¬ 
firming that including wave-function renormalization, in 
conjunction with a particular non-bare fermion-photon 
vertex, does not affect the theory in a significant way. 
The other important result is the existence of a possibly 
critical fermion flavour number of roughly 3.35, which 
is also consistent with some theoretical expectations and 
other numerical results. 

It is the first time that such a study includes the mo¬ 
mentum and energy dependence of the fermion and pho¬ 
ton self-energies and of wave-function renormalization in 
a gap equation with a non-bare vertex, and this allows 
a more reliable description of the behaviour of the the¬ 
ory. We estimate the numerical uncertainty for the values 
quoted, which comes mainly from the convergence crite¬ 
ria imposed, at about ±10%. The next step for future 
studies should be the inclusion of the imaginary parts 
of the self-energies in the equations and the relaxing of 
the approximation A(p) = B(p), which could in principle 


influence the results of this investigation. 
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